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Abstract. We present a theoretical description of the thermal fluctuations in a solid-supported stack of 
lipid bilayers, for the case of vanishing surface tension 7 = 0 and in the framework of continuous smectic 
elasticity. The model is successfully used to model the reflectivity profile of a thin (16 bilayers) DMPC 
sample under applied osmotic pressure and the diffuse scattering from a thick (800 bilayers) stack. We 
compare our model to previously existing theories. 


PACS. 61.10.Kw X-ray reflectometry (surfaces, interfaces, films) - 87.16.Dg Membranes, bilayers, and 
vesicles - 87.15.Ya Fluctuations 


1 Introduction 

Lipid bilayer systems have been extensively studied us¬ 
ing X-ray techniques, including small-angle scattering and 
reflectivity measurements 11 . Among the various experi¬ 
mental configurations employed, the case of solid-support¬ 
ed stacks of bilayers is very important from a practical 
point of view, as these samples are amenable to power¬ 
ful, interface-sensitive scattering techniques (specular and 
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nonspecular scattering, grazing incidence difraction) avoid¬ 
ing the ambiguities associated with powder averaging i2] 

El- 

Besides the fundamental physical problems they raise, 
such as the molecular forces responsible for the inter-bi¬ 
layer interaction potential (steric, electrostatic, van der 
Waals etc.); the influence of static defects; the loss of long- 
range order in low-dimensional systems and the role of the 
fluctuations, these systems can also serve as testing ground 
for the interaction of the cell membrane with membrane 
active molecules, such as antimicrobial peptides [4]. It is 
then of paramount importance to understand the differ- 
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ent parameters influencing the scattering signal in order to 
achieve a comprehension of the spectra before the (some¬ 
times subtle) effect of the included molecules can be con¬ 
fidently assessed. 

One of the most characteristic hallmarks of lamellar 
systems (exhibiting one-dimensional order) is the Landau- 
Peierls effect, whereby the long-range order is destroyed 
by thermal fluctuations. Ever since the seminal paper of 
Caille [5], this phenomenon has been studied in great de¬ 
tail, first in the bulk [M7ll8ll MTnilTrilT2] and subsequently 
in thin films mainly under the influence of 

extensive experimental studies on freely-suspended films 
of thermotropic smectics (see [T7J for a review). Recently, 
theoretical models were developed for the study of fluctu¬ 
ations in smectic films on solid substrates mm, taking 
into account both the surface tension at the free surface 
and the boundary conditions imposed by the substrate. 

In this paper we describe the thermal fluctuations in a 
solid-supported stack of lipid bilayers, concentrating on 
the (experimentally relevant) case of vanishing surface 
tension at the top of the stack. The influence of the sub¬ 
strate is only considered inasmuch as it limits the fluctu¬ 
ations of the bilayers, neglecting direct interactions (e. g. 
van der Waals, electrostatic) which can become important 
in the case of very thin films. 


along the normal direction. In contrast with (free-standing 
or supported) films of thermotropic smectics, for fully hy¬ 
drated lipid multilayers the surface tension at the free sur¬ 
face vanishes. This greatly simplifies the theoretical treat¬ 
ment of the fluctuations. We conclude this section by a 
discussion of the in-plane variation of the correlation func¬ 
tion. 

Section [3] starts with a discussion of the specular struc¬ 
ture factor S(q z ). We then consider the influence on S 
of the coverage rate, which can vary due to the prepara¬ 
tion technique or to partial dewetting upon hydration. We 
compare our model with experimental data on the spec¬ 
ular scattering and then with an estimate of the in-plane 
correlation function obtained from the diffuse scattering 
signal. 

2 Smectic elasticity 

We consider in the following a solid-supported sample with 
periodicity d , of thickness L = Nd and extending over a 
surface S in the plane of the layers. We take the origin of 
the z axis on the substrate, so that z = L gives the position 
of the free surface. The in-plane position is denoted by 
r_L = (x,y). 

2.1 Model and fluctuations 


The paper is structured as follows : in section [2] we 
discuss the positional fluctuations of the layers in a solid- 
supported stack in the framework of continuous smectic 
elasticity, taking however into account the discrete nature 
of the system by the limitation of the number of modes 


The simplest description of smectic elasticity is provided 
by the continuous smectic hamiltonian : 


F = \ [dr 


B 


( du(r ±1 z) 
V dz 


+ K ( A ± u(rj_,z)f 


+ 2 / drj " ' 


(1) 
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Following the treatment of Poniewierski and Holyst 
[ 20 ] , we shall decompose the deformation over independent 
modes; first, we take the Fourier transform of u(r±,z) in 
the plane of the bilayer : 


u{r±, z) = ~7=YL exp(-iqj.r±)w(q_L, a;), (2) 

V,b q_L 

The free energy m can now be written as the sum F = 

E qx F qm with ; 

1 


^q± — 


f dz 

B 

du{q±,z ) 

Jo 


dz 


+ k q± |w(q_L,^)r 


+ ^qi Kqj.,-01 2 . 


( 3 ) 


The boundary conditions for the tt(q_L, z) components are : 


ditions © : 


N 

u(q±,z) = ^2 Su n (qj_)f n (z) (7) 

n= 1 

where the summation goes from 1 to N, instead of oo, 
because only N components are required to describe the 
position of the N bilayers (this amounts to restricting the 
summation to the first Brillouin zone). Keep in mind, how¬ 
ever, that the index n denotes here a particular deforma¬ 
tion mode and not an individual bilayer. 

Finally, the free energy can be written as a sum of inde¬ 
pendent modes : F = \ J2q ± En=i A ( n > qjJ |^n(qj_)| 2 , 
where the ’’stiffness” A(n, qj_) associated to each mode 
depends on the elastic constants B , K and 7 . The equipar- 
tition theorem yields : 


iq±_u(q±. L) + B 


w(q_L, 0) = 0 
du(q±,z) 


dz 


= 0 . 


(4a) 

(4b) 


:=L 


The first condition © 1 ) simply states that the fluctuations 
go to zero at the substrate; the second one ©>) is neces¬ 
sary in order to write F q in © as a quadratic form [TdQ • 
Physically, it expresses the continuity across the interface 
of the a zz component of the stress tensor. 

The correlation function of the fluctuations can be de¬ 
fined as : 


z') = (u(r±, z)u{ 0, z')) , (5) 

where (•) denotes the ensemble average. From the Wiener- 
Khinchin theorem, its Fourier transform is : 


{Su n (q±)S 'Urn ( qjJ) 


= 6 


mn 


(| du n (q_L) | 2 ^ 


— 3mnk bF jA(ji , q_r), 


( 8 ) 


with S mn the Kronecker symbol. Plugging © in © , one 
has : 


C(q±,z,z') = ^rfn(z)f n {z')(\5u n (q 1 _)\ 2 'j . (9) 

n= 1 

We can now Fourier transform back to the real space do¬ 
main : 


C(r, z, z’) = F fn{z)f„{z') 

n =1 


( 10 ) 


q±dq±J 0 (qj_r) ^|<5u n (q_i_)| 2 ^ . 
where r = |rj_| is the in-plane distance (the _L symbol can 
safely be omitted). 



C(qj_, z, z r ) = (tt(qj_, ;;)u(-qj_, z’)) . (6) 2.2 The case of vanishing surface tension (7 = 0) 


We then expand u(q±, z) over the orthonormal set of har- Until now we have presented the general formalism; in 
monic functions f n (z), chosen to fulfill the boundary con- each case, we must find the orthonormal set of functions 
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f n (z), which are selected by the boundary conditions, and 
then determine the amplitude of each mode ^|<5« n (qj_)|"^ 
(or, equivalently, the stiffness A(n, qj_)). This is what we 
shall now do for the case when 7 = 0 , which is not only 
the simplest, but also the one relevant for systems of fully 
hydrated membranes m- 

The boundary conditions a are in this case : u(q±, 0 ) = 
= 0 , so the set of f n (z) is : 


0 and 




/ 2 . z 

/„(*) = Jjsm I 


(ID 


The amplitudes are given by : 

k B T 


l<5Wn(q-L)| 2 ) = 


B 

■B 


2n-l ny + ^ q 4 

4L 1 2 1 


2 L 

k B T A T 2 


( 12 ) 


B (2 n - 1) 2 tt 2 1 + (£ n q±) 4 


where £ 2 = 


2LX 


(2 n — 1 ) 7 r 


with A = \JK/B. We also define 


the dimensionless parameter 77 = — " B , first introduced 

2 B\d z 

by Caille [1 . 

Finally, from equation m we have : 

2 N A 


C(r,z,z')=v[-) J2 


7 r 7 ' ( 2 n — 1 ) 7 r 

7 n—1 v 7 


' 2 n —1 . (2n—l z'\ ( r 


(13) 


M(x) 



Fig. 1. The M(x) function 


where the function A4 can be expressed in terms of the 


Meijer G functioq 


~ 4 °°^ (^56 I 0 , 5 ’ 5 ’°) ( 14 ) 


A first observation is that the fluctuation amplitude 
C(0,z,z) = ^|u(0,i:)| 2 ^ is always finite ; the physical rea¬ 
son is that the presence of the substrate sets a lower 
boundary on the value of the wave vector in the z direc¬ 
tion : from equation m, q z > 7t/ 2L, thus forbidding the 
soft mode with q z —> 0 and suppressing the Landau-Peierls 
instability. Hence, there is no need for a lower cutoff in the 
integral in equation ra- 

it is immediately obvious that Ad(0) = 7 r/ 4 , so the 
correlation function for r± = 0 reduces to the simple for¬ 
mula : 

2 N 


c <o- 

x 7 n—1 


2n — 1 z \ . (2n — 1 z' 

sm | —-— 7 r— 1 sin 


(15) 


2 n L 


which we shall use in determining the specular scattering 
of the sample (subsection 13.11) . It is noteworthy that this 
function varies as C(0,z/d,z'/d) ~ ^(d/n) 2 for a fixed 
number of layers, so that any change in the thermody¬ 
namic parameter rj only results in a scale factor. This is 
very convenient for fitting the scattering spectrum of the 
system (see section [3]), since the time-consuming calcula¬ 
tion of the correlation matrix between the layers must only 
be performed once, and adjusting the 77 and d parameters 
only changes a prefactor. 

1 It can be represented in Mathematica by the sequence : 

iMeijerG [{{} , {}} , {{0, ±, A} , {0}} , (f) 4 ] 
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As an illustration, we present in Figure [2] the values of 
the correlation function in a stack of 100 bilayers; in (a), 
C{ 0, z, z ) represents the fluctuation amplitude for each bi¬ 
layer, while in (b) (7(0, z, L/2) represents the correlation 
of each bilayer with the one in the middle of the stack. 
Note the sigmoidal shape of the function in the first case 
and the very sharp peak in the second one. 




Fig. 2. The correlation function (scaled by r/(d/n) 2 ) in a 100 
layer film a) (7(0, z, z) b) (7(0, z, L/2). 


Although for any finite number of layers TV the fluc¬ 
tuation amplitude of the top layer (7(0, TVd, TVd) remains 
finite, it should diverge in the bulk limit TV —> oo. From 


formula m one obtains : 


< 7 ( 0 , Nd, Nd) 

V(d/ A ) 2 


N 


E 


i 

2n- 1 


1 

2 


7 + In 4 + ip 



(16) 


where 7 = 0.5772... is Euler’s constant and the digamma 
function ip(z) is the logarithmic derivative of the gamma 
function, given by : ip(z) = r'(z)/r{z). <7(0, TVd, TVd) 
diverges logarithmically, if (TV — 1/2) being indistinguish¬ 
able from In TV as soon as TV > 3. 

The same divergence spuriously appears even for a fi¬ 
nite L when the stack is treated as a completely continuous 
medium, without limiting the number of q z modes to TV 
(formally, this amounts to letting TV —> 00 and d —> 0 at 
fixed L = Nd). In this case, after summing over the q z 
modes, an artificial lower cutoff has to be introduced in 
the integral over q±, as in reference 114) . A comparison 
between our results and those obtained by this method is 
presented in the Appendix. 


2.3 In-plane variation 


We shall now consider the r variation of the correlation 
function. At this point, it is convenient to introduce the 
correlation of the height difference g(r,z, z') defined by : 

9 (r,z,z') = ((u(r ± ,z) - w(0,;?')) 2 \ 

X ’ (17) 

= (7(0, 0 , z) + <7(0, z', z') - 2<7(r, 0 , z'). 

which has the advantage of remaining finite (for finite val¬ 
ues of r) even in unbound systems. We shall further write 
the argument of the A4 function in equation m as : 

I f 277- 117T 7* 

Y--, where £ = \f\d emphasizing that, for a 
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given number of layers, the r variable in C{r,z,z') and 
g(r, z, z') scales with the correlation length £. The phys¬ 
ical significance of this quantity can be seen as follows 
[22] : for distances less than £, the layers fluctuate inde¬ 
pendently, while for distances greater than £ the fluctua¬ 
tions are coherent from layer to layer. For an unbounded 
medium (translation invariance along z) it was shown |15| 
that, for z = z' the g(r) function behaves as 


77 w He 


,\ 2 i , \ 2 

a \ 1 r 


1 — 7 +In 2 — In ( — 


(18) 


if r < £ and as 


r 




In ( - 


(19) 



Fig. 3. The height difference self-correlation function 
g(r/£, z, z) scaled by g(d/n) 2 in an 800 layer film, for z = lOOd 
(triangles) and z = 800d (circles). At small values of r/£, both 
curves are well described by the asymptotic form m deter¬ 
mined by Lei [15] (solid line), with no adjustable parameters. 
Also shown is the average value g avs {r) (dotted line). 


if r > £. In the limit r —> oo, g(r) diverges logarithmically, 
as it is well-known from the continuum theory [5). 

We present in Figure[3]the height difference self-corre¬ 
lation function g(r/£, z, z), computed from equations (1171) 
and m for a stack of 800 bilayers. 

In the low r limit we find that, for all values of z , g(r) is 
very well described by the asymptotic form (H51) . For high 
values of r/£, the divergence (THU) is replaced by saturation 
to a value of g( i —> oo) = 2(7(0). 

As we shall see in subsection 13.51 the self-correlation 
function g(r ) can be extracted from the diffuse scattering 
signal; however, when g also depends on 2 this signal will 
be an average over the positions in the stack. For com¬ 
parison with the experimental data we calculate ga, vg (r) 
(shown in Figure [3] as dotted line) as an average over 
g(r, z,z) for eight different values of 2, equally spaced from 
100d to 800d. 


3 Scattering 

As discussed in the Introduction, notable differences ap¬ 
pear between the scattering signal from bulk smectic phases 
and from solid-supported films. Even if the latter contain 
thousands of bilayers, the clear separation between spec¬ 
ular and diffuse scattering shows that the Landau-Peierls 
effect is suppressed [231 . 

3.1 Specular scattering 

In specular reflectivity studies, the incidence and reflection 
angles are equal : = a/, corresponding to q± = 0. As 

we shall see below, in this case only the correlation at 
zero in-plane distance (given by formula (1151) ') is relevant. 
For simplicity and to emphasize the discrete nature of the 
stack we shall use the notation (7(0, nd, md) = <7(n, to). 
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The specular scattering factor of the lamellar stack 
(without taking into account the substrate) S(q z , qj_ = 0 ) 

can be written as mm ■ 

N 

S(q z ) = ^ cos [q z d(m-n)]e-^^ c{rn ’ m)+G ^' n)) ( 20 ) 

m,n =1 

The spectra are corrected for the diffuse scattering by 
subtraction of an offset scan [ 2 ]. 



approaching the free surface. The structure factor then 
reads : 


N 

S{q z )= ^2 f{n)f(m) ■ cos[q z d(m-n)} ( 21 ) 

m,n—1 

q 2 

•exp C(m,m) + C(n,n )) 


For the coverage function we chose the convenient an¬ 
alytical form f(n) = [l — ( 77 )°] , where a is an empiri¬ 
cal parameter controlling the degree of coverage. This is 
a convenient method, but not a very precise one, insofar 
as the fluctuation spectrum is still calculated for a fixed 
number of layers, N. In the limit where the size of the 
domains with a given thickness is larger than the X-ray 
coherence length, an alternative approach would be to in¬ 
coherently (no cross-terms between domains) average over 
a distribution depending on the total layer number P(N). 


Fig. 4. Specular structure factor for a stack of 1000 bilayers, 
with d — 62 A, for different values of the Caille parameter 
(from top to bottom, tj = 0, 0.1, 0.2, and 0.3). Curves vertically 
shifted for clarity. 


3.2 Accounting for partial dewetting 

Fully-hydrated lipid multilayers differ from freely-suspend¬ 
ed smectic films in that their thickness is not constant; 
partial dewetting leads to a variation in thickness over the 
area of the sample m- Notably, this leads to an obliter¬ 
ation of the Kiessig fringes, due to interference between 
the external interfaces of the system. We take this ef¬ 
fect into account by introducing a coverage rate for each 
layer, f(n), which is 1 at the substrate and decreases on 


3.3 Describing the reflectivity profile 

In order to model the reflectivity profile of our system, 
besides the structure factor of the stack one needs the form 
factor of the bilayers. Furthermore, the presence of the 
substrate must be taken into account. This is done using 
a semi-kinematic approximation, where the reflectivity of 
a rough interface is expressed by the master-formula of 
reflectivity m ■■ 


R{q z ) = Rpiqz) ■ 


-f 

P 12 Jo 


dp(z) 

dz 


e iq * z d z 


( 22 ) 


where I?f denotes the Fresnel reflectivity of the sharp in¬ 
terface and p(z ) is the intrinsic electron density profile 
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whereas p \2 is the total step in electron density between 
the two adjoining media (silicon and water, in our case). 

For the electronic density profile of the bilayer we use 
a parameterization in terms of Fourier components [2], 
which has the advantage of describing the smooth pro¬ 
file of lipid bilayers using only a few coefficients. For one 
bilayer we have : 



where Ph\(z) is defined between —d/2 and d/2, N comp is 
the total number of Fourier components, pk is the ampli¬ 
tude of the fc-th component and Vk the associated (com¬ 
plex) phase factor, which in our case can be shown to 
reduce to ±1 only, due to the mirror symmetry of the 
bilayers. The bilayer form factor is given by the Fourier 
transform of pb\(z). 

The total density profile is given by the density pro¬ 
files of the N bilayers (weighted by the coverage factors / 
described in subsection to which is added the profile 
of the substrate, described by an error function of width 
the rms roughness of the Si wafer (this quantity can be in¬ 
dependently determined from the reflectivity of the blank 
wafer and is typically worth 8 — 10 A for all our measure¬ 
ments). 

Effects related to the finite instrumental resolution and 
sample absorption are completely negligible in thin sam¬ 
ples and were not implemented. 


3.4 Comparison with experimental data 

Lipid films partially hydrated in humid atmosphere ex¬ 
hibit many well-defined Bragg peaks, indicative of very 
low fluctuation amplitudes, so they are not a good test¬ 
ing ground for our model. On the other hand, samples hy¬ 
drated in excess pure water sometimes have too few Bragg 
peaks. We thus chose to test our model against spectra ob¬ 
tained on samples that are in excess solvent, but under an 
osmotic pressure imposed by PEG (polyethylene glycol) 
solutions [2Sj. Fig. [5] shows the X-ray reflectivity of 16 
l,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bi¬ 
layers on a silicon substrate measured in an aqueous PEG 
solution of 3.6 wt. % concentration. The temperature was 
maintained at 40 °C, at which the lipids are in the fluid L a 
phase. The continuous line shows a full g-range fit to the 
data using the structure factor (ED- From the position of 
the Bragg peaks one can determine a periodicity of 59.5 A. 
The suppression of the higher-order Bragg peaks clearly 
indicates the influence of the fluctuations. The fit yielded 
the parameters: p = 0.065 for the Caille parameter and 
a = 1.7 for the coverage exponent in 1211 The correspond¬ 
ing fluctuation amplitudes cr 2 (n) = C(n,n) and coverage 
function f(n) are also shown. Further details on the re¬ 
sults of these measurements will be presented elsewhere 

ESI- 

Experimental details : The curve in Fig. [5] was mea¬ 
sured at the bending magnet beamline D4 at the HASY- 
LAB/Desy in Hamburg, Germany. After the beam passed 
a Rh-Mirror to reduce high energy flux, a photon en¬ 
ergy of 19.92 keV was chosen using a single Si(lll) crys- 
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Fig. 5. X-ray reflectivity curve of 16 DMPC bilayers on a Si substrate measured under osmotic pressure (symbols) fitted with 
a full g-range semi-kinematic model [2] using the structure factor (12111 (solid line). Inset : (top) fluctuation amplitude cr 2 {n) for 
a Caille parameter r/ = 0.065 obtained from the fit and (bottom) coverage function /(n) with a = 1.7 deduced from the fit. 


tal monochromator. The beam was collimated with sev¬ 
eral motorized vertical and horizontal slits. The reflected 
intensity was measured with a fast scintillation counter 
(Cyberstar, Oxford). 

Sample preparation: The lipid 1,2-dimyristoyl-sn- 
glycero-3-phosphocholine (DMPC) was bought from A- 
vanti and used without further purification. 16 bilayers 
were prepared on a commercial silicon substrate using the 
spin-coating method as described in 27]. The substrate 
was cut to a size of 15 x 25 nun 2 and carefully cleaned 
in an ultrasonic bath with methanol for ten minutes and 
subsequently rinsed with methanol and ultrapure water 
(Milli-Q, Millipore). The lipid was dissolved in chloroform 
at a concentration of 10 mg/ml. An amount of 100 fjl of 


the solution was pipetted onto the cleaned and dried sub¬ 
strate which was then accelerated to 3000 rpm using a 
spin-coater. After rotation for 30 seconds, the sample was 
exposed to high vacuum over night in order to remove all 
remaining traces of solvent. The sample was refrigerated 
until the measurement. For the measurement the sample 
was mounted in a stainless steel chamber with kapton win¬ 
dows which can be filled and in situ flushed with fluids 
such as water or polymer solutions. PEG of molar weight 
20.000 was bought from Fluka and used without further 
purification. A concentration of 3.6 wt. % corresponds to 
an osmotic pressure of about 10 4 Pa [25] , 
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3.5 Diffuse scattering 

It is well-known that the diffuse scattering signal from 
a surface can be expressed in terms of the height dif¬ 
ference self-correlation function g{r) of the surface [ 25 ,. 
For a multilayer system, one must also take into account 
the cross-correlation function g(r,i,j ) [30] and the calcu¬ 
lations become much more complicated. However, it can 
be shown that by integrating the diffuse signal in q z over 
a Brillouin zone, the cross-correlation terms i ^ j cancel 
and one is left with a curve corresponding to a transform 
of an average self-correlation function m- 

We compare our model to recent experimental data 
obtained on stacks of 800 fully-hydrated DMPC bilayers, 
where the experimental self-correlation function g ex p ( 7 ') 
was obtained by inverting the integrated diffuse scattering 
[32]. In Figure [G] we present g ex p (r) (open dots), as well 
as the ga, vg (r) function defined in subsection !2.3l with two 
different sets of fitting parameters r) and if. 

The theoretical curve can be scaled to the experimen¬ 
tal data, either by adjusting both parameters, yielding 
Tj = 0.071 for the Caille parameter and £ = \f\d = 20 A 
(solid line in Fig. [ 6 ]) or with a fixed parameter £ = 40 A 
(from the analysis of the half width of the diffuse Bragg 
sheets [32]), and an open parameter ry which is adjusted 
to 77 ~ 0.09. Note that the present expression for g{r) 
can give better account of the experimental data than the 
bulk correlation function used before 0332 ]. In particu¬ 
lar, it reproduces the observed saturation regime at high 
r, while the bulk correlation function diverges logarithmi¬ 
cally. However, systematic discrepancies between data and 



Fig. 6. Experimental values for the self-correlation function 
9exp (r) (open dots) and < 7 avg (r) for two sets of fitting param¬ 
eters : g = 0.071 and £ = 20 A (solid line) and g = 0.09 and 
£ = 40 A (dashed line). Inset : the same graph in log-log scal¬ 
ing. 

theory should not be overlooked. We note that the func¬ 
tional form of g[r) for small r <C £ at is not well captured, 
in particular it can not explain the exponent g(r) oc r 0 ' 7 
of the algebraic regime in the data. The possible reasons 
for this discrepancy can be manifold: contributions of non¬ 
bending modes to the diffuse scattering, a length scale de¬ 
pendent bending rigidity k, residual tension in the bilayers 
due to edge effects, or nonlinearities in the Hamiltonian. 
This needs to be investigated in future studies. Finally, a 
more detailed treatment of the 2 -dependence is needed, 
since the experimental analysis determines a correlation 
function which is averaged over the scattering volume. 

4 Conclusion 

The presence of a substrate dramatically changes the ther¬ 
mal fluctuations in lipid multilayers. Most noticeably, the 
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Landau-Peierls instability is suppressed. We present a the¬ 
oretical model taking into account this feature and show 
that it describes very well the experimental reflectivity 
data. Reasonable agreement is also obtained for the dif¬ 
fuse scattering. Our result is a first step towards a unitary 
interpretation of the global (specular and diffuse) scatter¬ 
ing signal of solid-supported lipid multilayers. 

Appendix. Comparison with previous results 

In order to assess the validity of our method, we compared 
our results with those obtained by Romanov and Ul’yanov 
by a rigorous treatment of the discrete model [T5]. How¬ 
ever, they only show data for smectic films with a rather 
high surface tension (7 = 30mN/m). Thus, we chose their 
thickest film (21 layers) and only compared the values of 
C for the bottom half (close to the substrate), where the 
surface tension plays a lesser role. As shown in Figure [3 
the agreement is quite good. 

We also compared our results to those obtained by 
Shalaginov and Romanov )T4] , who also used a continuous 
model, but without restricting the number of modes. For 
a solid substrate and a vanishing surface tension at the 
top, their equations (18a) and ( 22 ) lead to : 

k B T dx 

G (m,n) = - ; - / -—— 

8itVBK Jo x cosh(a;JV) 

[sinh[a;(n + m — N )] + sinh(xlV) cosh[x(ra — m)\ 

— cosh (xN) sinh(x |n — m |)] 

In Figure [5] we present our results for C(n,n) in a 
film with 61 layers, with d = 30 A and rj = 0.14 (solid 
dots) as well as the values obtained using the formula (HHl) . 



Fig. 7. Comparison between the results of the present method, 
for 7 = 0 (in full symbols and solid line) and the values ob¬ 
tained by Romanov and Ul’yanov m for 7 = 30 mN/m (open 
symbols, dotted line) for the correlation function in a film com¬ 
prising 21 layers (with d = 30 A and rj = 0.14). Values only 
shown for the bottom half of the stack. 


for different values of the cutoff parameter a 0 (from top 
to bottom, 1.5, 4, 10, 30, 85, and 200 A). Clearly, the 
fluctuation amplitude is very sensitive to the value of a 0 , 
with an approximately logarithmic dependence, and gives 
the same results as our method for cio = 85 A. We remind 
that neither in our method, nor in the one of Romanov 
and Ul’yanov is there any need for a cutoff, due to mode 
number limitation (see the discussion in subsection 12 .211 . 

D. C. has been supported by a Marie Curie Fellowship of 
the European Community programme Improving the Human 
Research Potential under contract number HPMF-CT-2002- 
01903. 
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Fig. 8. Comparison between the results of the present method 
(solid dots) and the values obtained by the method of Shalagi- 
nov and Romanov [14] (lines), for a film comprising 61 layers. 
From top to bottom, the values of the cutoff parameter ao in 
equation (l2dl) are : 1.5, 4, 10, 30, 85, and 200 A. 
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